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We present a detailed description of high-energy neutrino and antineutrino inelas- 
tic inclusive scattering off nuclei in terms of nuclear structure functions. In our 
approach we take into account a QCD description of the nucleon structure func- 
tions as well as a number of basic nuclear effects including nuclear shadowing, Fermi 
motion and binding, nuclear pion excess and off-shell correction to bound nucleon 
structure functions. These effects prove to be important in the studies of charged- 
lepton deep-inelastic scattering. We discuss similarities and dissimilarities in the 
calculation of nuclear effects for charged-lepton and neutrino scattering caused by 
nonconserved axial current in neutrino scattering. We examine the Adler and the 
Gross-Llewellyn-Smith sum rules for nuclear structure functions and find a remark- 
able cancellation between nuclear shadowing and off-shell corrections in these sum 
rules. We present calculations of differential cross sections for inclusive neutrino and 
antineutrino scattering in comparison with recent data on different target materials. 

PACS numbers: 13.60.Hb, 25.30.Pt, 24.85. +p 

I. INTRODUCTION 

The interest in the study of neutrino interactions with nuclei at intermediate and high 
energy has grown considerably in recent years (see, e.g., [lj). New experimental results on 
neutrino cross-sections are also obtained for different nuclear targets 0, 0, 0| ■ The presence 
of an axial-vector component in the weak current and the quark flavor selection distinguish 
neutrino interactions from charged-lepton (CL) and hadron collisions. This feature makes 
the neutrino (is) and antineutrino (is) data a unique source of information on the nucleon and 
nuclear structure. The understanding of neutrino propagation in matter is also important 
for astrophysics, cosmology and even geology applications. 

The use of heavy nuclear targets in neutrino experiments is motivated by the need to 
collect a statistically significant number of interactions. For this reason, the understanding 
of nuclear effects is of primary importance for a correct interpretation of experimental results 
and for the evaluation of the corresponding uncertainties. An example of such synergy is 
provided by the precision electroweak measurements from neutrino interactions. The role of 
nuclear corrections to neutrino structure functions has recently been emphasized jl| in the 
context of the anomalous measurement of the weak mixing angle (sin 2 9w) reported by the 
NuTeV collaboration 0|. 

An accurate account of nuclear effects is not only important in the determination of 
electroweak parameters in neutrino scattering experiments, but also for the understanding 
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of neutrino masses and mixing. The next generation experiments would imply precision 
measurements to disentangle small oscillation signals from neutrino and antineutrino inter- 
actions on nuclei. This in turn would require us to improve our knowledge of i^(i>)-nucleus 
cross-sections in order to reduce systematic uncertainties. 

The availability of high- intensity neutrino beams from the recent NuMI {?} and JPARC [sj 
facilities offers new opportunities for detailed studies of cross sections and nuclear effects 
in neutrino interactions. One such example is MINERz/A {9J, a dedicated experiment at 
Fermilab which will collect data in few years. On a longer time scale, the construction of a 
neutrino factory [Io[ would then allow a further step forward, finally reaching the ultimate 
precision of the neutrino probe. 

In this work we present a calculation of inelastic neutrino-nucleus structure functions 
(SF) and differential cross sections in a wide kinematical range of the Bjorken variable x and 
momentum transfer squared Q 2 . Our goal is to develop a quantitative model incorporating 
existing data, which would be useful in interpreting v(v) experiments. 

Scattering experiments with charged leptons show that even in the deep-inelastic scatter- 
ing (DIS) region, for which the energy and momentum transfer are larger than the nucleon 
mass M, nuclear SF differ from a simple sum over bound nucleons and a significant nuclear 
effect is present even in the scaling regime (see, e.g., [H[). As evidenced by numerous studies 
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12J, ll3| and also [14(), the lepton-nucleus DIS is characterized by different mechanisms 



in different kinematical regions of the Bjorken variable x. In a recent paper 14| we studied 
the charged-lepton DIS off nuclear targets and developed a quantitative model for nuclear 
structure functions, taking into account major nuclear effects including nuclear shadowing, 
Fermi motion and binding, nuclear pion excess and off-shell correction to bound nucleon 
structure functions. This approach showed a very good agreement with data on the EMC 
effect for light and heavy nuclei for both the x and Q 2 dependence. In the present paper we 
extend this approach to neutrino-nucleus inelastic scattering. 

Although nuclear scattering mechanisms for high-energy charged leptons and neutrinos 
are similar in many respects (especially in the impulse approximation), there are important 
differences due to the presence of the axial- vector current in neutrino interactions. The inter- 
ference between the vector and the axial-vector currents introduces C-odd terms in neutrino 
cross sections, which are described by SF F 3 . In the calculation of nuclear corrections, we 
separate the contributions to different SF according to their C-parity. We examine in detail 
the dependence of nuclear effects on C-parity and show that it is particulary important in 
the nuclear shadowing region. Results are then used to predict nuclear structure functions 
for neutrino and antineutrino interactions, respectively. 

In contrast to the electromagnetic current, the axial current is not conserved, and for 
this reason the neutrino structure functions in the region of low Q 2 and small x are rather 
different from those of electroproduction. In particular, the neutrino longitudinal SF F^ 
(as well as F 2 ) does not vanish at low Q 2 . Its value is determined by the divergence of the 
axial- vector current which is linked to the virtual pion cross-section (Adler relation [TEI). We 
discuss the derivation of this relationship for structure functions and examine corrections to 
that. We also address the phenomeno logical implications of this relation in the description 
of low-Q 2 and low-x neutrino differential cross sections. The nuclear shadowing effect on the 
virtual pion cross section results in a significant suppression of the corresponding nuclear 
structure functions. 

The paper is organized as follows. In Sec. [Til we introduce the general formalism of neu- 
trino inclusive inelastic scattering, define the structure functions and present the differential 
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cross sections for charged- current (CC) and neutral-current (NC) neutrino scattering. In 
Sec. Ill CI we summarize our model of the nucleon SF in the region of high Q 2 ; in Sec. IIIDI wc 
deal with the low-Q 2 region and discuss the separation of the axial current contribution from 
the vector component. In Sec. II III we discuss the nuclear effects in neutrino DIS. Particular 
attention is paid to the calculation of neutrino C-even and C-odd structure functions. In 
Sec. IIII Dl and IIIIEI we examine the Adler and the Gross-Llewellyn-Smith sum rules for 
nuclear targets. In Sec. [IV] we present the results of numerical analysis of SF and differential 
cross sections for a number of nuclear targets and compare our predictions with the recent 
data UMM- 



II. GENERAL FORMALISM OF NEUTRINO INCLUSIVE SCATTERING 

To leading order in the weak coupling constant, neutrino scattering is described by the 
standard one-boson exchange process. Neutrino interactions can be mediated by the charged 
W ± boson (charged curent or CC) or by the neutral Z boson (neutral curent or NC). In the 
Standard Model (SM) the leptonic and hadronic charged currents can be written as 

j~ = eO L x v e + fiO^, (la) 
J- = d'O^u + sO^c, (lb) 

where 0\ = (1 — 7s)7a, and d' = d cos 9c + s sin 9c, s' = —d sin 9c + s cos 9c are the super- 
positions of d- and s-quark states and 9c is a Cabibbo angle (sin 2 9c ~ 0.05). For simplicity, 
we neglect contributions from the third generation of quarks and leptons. The superscript 
L in Eqs.Q indicates that only left doublets participate in the CC weak interaction. 
The neutral current can be written as 

Ja = E^-^7 5 )7a^, (2a) 

t 

J°x = T,«(9 q v -9hshxq, (2b) 

where the sum is taken over all types of leptons (£ = e, u e , /i, v^) and quarks (q = u, d, c, s). 
The vector and axial charges of quarks and leptons in the SM in terms of the weak mixing 



angle 9w can be found in, e.g., [161 . 



A. Charged-curent neutrino scattering 

We first consider the CC (anti)neutrino inelastic inclusive scattering. In inclusive scatter- 
ing, the final hadronic state is not detected and the differential cross section can be written 
as (see, e.g., 

_ 2^/(47r)3 d*k> 
cc " (l + QVM 2 w Y L ^ X J^k)E" (3) 

where Gf is the Fermi weak coupling constant; Mw is the IV-boson mass; p is the target four- 
momentum; k = (E, k) and k' = (E', k') are four-momenta of the incoming and outgoing 
lepton, respectively; q = k — k' is four-momentum transfer and Q 2 = —q 2 . The tensors L^x 
and W^x describe the interaction of the W boson with leptons and hadrons, respectively. If 
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the polarization of the outgoing lepton is not detected, then the leptonic tensors for neutrino 
and antineutrino scattering are 

L ( X ] = § (k,k' x + k x k^ - k ■ k'g, x ±ie^(k, k')) , (4) 

where the sign +(— ) corresponds to the neutrino (antineutrino). We use the contracted 
notation e^a, b) = e^xa^b^ . The hadronic tensor W^x is the sum of CC matrix elements 
over all possible final hadronic states. For the neutrino we have 

w$>(p,q) ^^j2( 2 ^ 5 (p+9-Pf)wfmf)(f\jm\p), (5) 



where the current in the SM is given by Eq. fllbj) . The antineutrino tensor corresponds to the 
exchange of the W~ boson and is given by a similar equation with the current J x replaced 

by J \ ■ 

For the scattering from an unpolarized target, there are five independent structure func- 



tions in the hadronic tensor §5§ for either the neutrino or antineutrino (see, e.g., [16|) 



w^xip, q) = - g^xj f x + (p^ - q/~r) (px - 0a^t-J 
. . / n ^3 q^qx „ . q^px + qxp^ „ 

+ «£a*a(p,5)o + "7^ F 4 + — -F s . 

2p ■ q Q z p ■ q 



p-q 



(6) 



We use the normalization of states (p\p') = 2E p (2ir) 3 6(p — p') for both the bosons and 
fermions. In this normalization the hadronic tensor (JSJ) and the structure functions are 
dimensionless. The terms with F\ and F 2 are similar to those in charged-lepton scattering. 
They originate from vector- vector and axial-axial correlations in Eq.(jSJ). The antisymmetric 
term (F$) describes parity- violating vector-axial and axial- vector transitions. The terms 
and F5 are present because of nonconservation of the axial current. 

The neutrino and antineutrino structure functions are apparently different. We will also 
consider the sum and the difference of neutrino and antineutrino structure functions, which 
have definite C-parity 

F u±, = F u ± (?) 

Contracting the leptonic and hadronic tensors, we obtain the explicit form of the differen- 
tial cross section in terms of the structure functions. Using the usual variables x = Q 2 / (2p-q) 
and y = p ■ q/p ■ k we have 

^cc _ G 2 F (p-k)/n {v>v) 



dxdy (l + Q'/M^y^ 
The kinematical factors read as follows 



TYiF?*. (8) 





m \ 


Q 2 V 






\ 2 y 2 Q' 2 


2Q 2 , 


) 4Q 2 



y = y 2 ^ ( 1 - ^ ) > M 

P 2 Q 2 
(p ■ qf 



Y = (1-^1 -^(1 + 7^1. 
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yV 2 

2Q 



^3 = ±^(l-SU, (9c) 



/2 

m 



where the sign +(— ) refers to neutrino (antineutrino) scattering, m! is the mass of the out- 
going charged lepton and Q' 2 = Q 2 + m' 2 . We keep the lepton mass terms for completeness. 
Although these terms are tiny in v e and scattering, they are not negligible in v T scattering. 
The contributions from the structure functions F 4 and F 5 to the neutrino production cross 
section are suppessed at high energy by a small ratio m' 2 /p ■ k. 

The structure functions can be related to the virtual boson helicity cross sections by 
projecting Eq.([6]) onto the states with definite polarizations of an intermediate vector boson 
[for simplicity, we suppress explicit notation for (anti)neutrino hadronic tensor] 

W h = e^W^el (10) 

where is the polarization vector of the virtual boson in the helicity state h (h = ±1 
corresponds to two transverse polarizations, and h = to the longitudinal polarization). 
The relation between the helicity structure functions Wh and -^1,2,3 can be derived from 
Eq. ([6]) using an explicit form of polarization vectors and the orthogonality and normalization 
conditions for polarization vectors (see, e.g., [Ill, Oil). We have 



W± 1 = F 1 ±jF 3 , (11a) 
W = 1 2 F 2 /{2x)-F 1 . (lib) 

One concludes from Eqs. fllll) that Fi is the transverse structure function averaged over trans- 
verse polarizations and F 3 is determined by the right-left asymmetry in helicity structure 
functions W+\ - W-\. Below, we also use F T = 2xF t , F L = 2xW and R = F L /F T . 



B. Neutral-current neutrino scattering 

The NC hadronic tensor is given by Eq.fjolh in which the charged current must be re- 
placed by the neutral current. The leptonic tensor in this case is given by Eq.fll]) with the 
overall factor 1/2, reflecting the fact that the neutrino can only be left-polarized (and the 
antineutrino, right-polarized). As a result, the NC cross section has the overall factor 1/2 
compared to Eq.©, 

d 2 4? G 2 F ( P -k)/(2n) " 

dxdy (1 + Q 2 /M 2 ) 2 1 1 ' K ' 

where Mz is the mass of the Z boson and the kinematical factors Yi are obtained from 
Eqs.flS]) by taking the limit m' — > 0. Therefore, the structure functions F4 and F5 do not 
contribute to the NC cross section. We also comment that neutrino and antineutrino NC 
structure functions of the type i are identical, because the neutrino and antineutrino in NC 
scattering couple to the same hadronic current. 
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C. Neutrino structure functions at high Q 2 

At high Q 2 , above the resonance region, the structure functions are usually described in 
QCD as a series in Q~ 2 on the basis of the operator product expansion of the correlator 
of the weak currents (twist expansion) [13]. The leading twist (LT) terms correspond to 
quasifree scattering off quarks corrected for gluon radiation effects. In this approximation 
the structure functions are given in terms of parton distributions or PDFs (see, e.g., [l6|, If 



The PDFs are universal high- momentum transfer characteristics of the target. The PDFs 
have a nonperturbative origin and their dependence on the Bjorken variable x can not be 
calculated in perturbative QCD. However, the dependence of PDFs on momentum transfer 
is governed by the perturbation theory in strong coupling constant a s at the scale Q 2 
(DGLAP evolution equations (l9jj) . Currently, the coefficient and the splitting functions, 
which determine the QCD evolution, are known up to two-loop approximation (NNLO 
approximation) 0, 2~l| . 



Along with the LT terms, the higher order Q~ 2 terms have to be taken into account 
(higher twists or HT). The HT terms can be of two different sources: (i) the kinematical 
target mass terms which can be approximately resummed and absorbed into the LT terms 



according to [22j, |23[ and (ii) the dynamical HT terms which are related to quark-gluon 
correlation effects. In summary, the structure function of type i can be written as a power 
series 

Ffa Q 2 ) = if MC (x, Q 2 ) + H^/Q 2 + ■ ■ ■ , (13) 

where F™ c is the LT structure function corrected for the target mass effects and are 
the functions of x and Q 2 describing the strength of the HT terms of twist t. In this paper 
the target mass corrections are computed using the approach of Ref . [22| . In order to keep 
the correct threshold behaviour of structure functions for x 1, we expand F™ c in a 
power series in Q~ 2 and keep the terms to the order Q~ 2 

In the numerical analysis the LT structure functions are computed using both the PDFs 
and the coefficient functions to the NNLO approximation. We use the PDFs and HT of 
Ref. 24j, obtained from a new fits optimized at low Q 2 and including additional data. It 



should be noted that unlike the PDFs the HT terms are not universal and generally depend 
on the type of the structure function i as well as on the type of interaction. In order to 
evaluate HT in neutrino DIS, we assume that the ratio Hi/ F^ 1 - is similar for CL and neutrino 
structure functions as a working hypothesis. That allows us to use the phenomenological 
HT terms extracted from the fits to CL DIS data. We also assume no HT terms for F3 
in the calculations of nuclear structure functions and differential cross sections presented in 
Sec. El 

In the calculation of neutrino cross sections for the structure function F 5 we use the 



Albright-Jarlskog relation 2xF 5 = F 2 (27| . Recently, it was argued that this relation survives 



QCD higher order corrections and the target mass corrections [28J. For F 4 we will use the 



relation F 4 = F 2 /(2x) — Fx which replaces the Albright-Jarlskog relation F 4 = for massless 



quarks if the target mass corrections are taken into account [28 



D. Neutrino structure functions at low Q 2 



We note that even for high-energy scattering in the fixed-target experiments the events 
with small Bjorken x typically have the values of Q 2 below lGeV 2 . In this Section we discuss 
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the vector- and the axial-current contributions to neutrino structure functions at low Q 2 and 
outline the model which will be used in our studies of neutrino differential cross sections. 
In this context it is useful to review the requirements which follow from conservation of the 
vector current and partial conservation of the axial current. It follows from Eqs. (1101) and 
(J5]) that the structure functions Ft and Fi can be written as 

F t ,l = -QV,l, (14) 

7T 

where <Tt(l) refers to the total cross section induced by the transverse (longitudinal) com- 
ponent of the weak current J\, 

^ = h E( 27r ) 4 ^ + i-Pf) I (p\ e T, L ■ J(o) I/) | 2 , (15) 
3 f 

where j = ((p ■ q) 2 — p 2 q 2 ) 1 ^ 2 is the invarint flux of the virtual boson with four-momentum q 
on the target with four-momentum p (for the nucleon at rest j = M\q\, with M the nucleon 
mass) and the sum is taken over all possible final states. The transverse cross section does 
not vanish in the limit Q 2 — > 0; therefore Ft vanishes as Q 2 in this limit. This holds for 
both the vector- and the axial-current contributions. 

In the longitudinal channel the low-Q behavior of the vector- and axial-current parts are 
different. We first consider the contribution from the vector current. The conservation of 
the vector current (CVC) suggests q^W^ = for the vector-current part of the hadronic 
tensor. From this condition we conclude that F^ c vanishes faster than F^P at low Q 2 and 
pvc^pvc ^ Q 2 /q 2 , similar to the charged-lepton case. 

In contrast to the vector current, the axial current is not conserved and, according to the 
hypothesis of partial conservation of the axial current (PCAC), its divergence is proportional 
to the pion field flEl ]. 

dA ± = Uml^, (16) 

where m 7T is the pion mass and f n = 0.93m,,- is the pion decay constant and (p^ is the pion 
field in the corresponding charge state. The PCAC relation can be used to compute the 
axial current contribution to the longitudinal structure function F£ c in the region of low 
momentum transfer. In this context it is useful to explicitly separate the pion contribution 
to the axial current, 

A, = -f n d^ + 4, (17) 

where the operator A' describes the contributions from heavy hadron states, i.e. the contri- 
butions from axial-vector meson states, pn continuum, etc. It is important to notice that 



the pion derivative term does not contribute to the structure functions [29J . This is because 
the matrix elements of d^ip are proportional to the momentum transfer q^. These terms 
give a vanishing contribution when contracted with the boson polarization vectors because 
of the orthogonality condition • q = for any helicity state h. Therefore, the axial- vector 
operator A^ can be replaced by A'^ everywhere in Eqs. §5$). The PCAC relation for A' follows 
from Eq. ffTB"]) : 

dA!± = fj±, (18) 

where j^r = (d 2 + m 2 ^^ is the pion current. In order to calculate Fl we separate the 
contribution of the axial current to the hadronic tensor (J3J), contract it with the polarization 
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vector eo and replace A with A'. 1 Using Eq. flIB"]) we have 

F£ c = 7 3 Ff AC + 27 3 ^ 2/ "^ + F£ c , (19) 

7r|q| 

where F^ CAC = ] 2 o^j-n and = cr n (s, Q 2 ) is the total cross section for scattering of a virtual 
pion with four- momentum q and the center-of-mass energy squared s = (p + q) 2 

3 f 



The flux j is defined in Eq. (fT5l) . Note also that ir + cross section corresponds to the neu- 
trino, and ir~ to the antineutrino. The quantity a' n describes the interference between the 
divergence of the axial current and the operator A' z in the hadronic tensor 

< = rS (27r)45(p+g " P/,)Im ( (p| ^ (0)l/,)(/,| ^ (0)tb) ) • (21) 

3 r 

The last term in Eq. (I19|) is similar to and vanishes as Q 4 . 

It follows from the relation F 2 = (Fl + Ft)/^ 2 and f)19p that the structure function F 2 
at low Q 2 is dominated by the i^ CAC term. The PCAC contribution to the neutrino cross 
sections at small values of Q 2 was experimentally tested in [3oT ]. On the other hand, the 
PCAC term should vanish at high Q 2 . In order to interpolate between low and high Q 2 , we 
introduce a form factor f PCA c(Q 2 ) 

f pcac = f>JL fpcAcXQ 2y (22) 

7T 

In the numerical analysis we assume a dipole form f PC Ac(Q 2 ) = (1 + Q 2 / M 2 CAC )~ 2 with M PCAC 
the mass scale controlling the PCAC mechanism. Since the pion pole does not contribute 
to the structure functions, the scale M PCAC is not determined by the pion mass, but rather 
it is related to higher mass states like a\ meson, pir continuum etc. In the present paper we 
fix the numerical value of M PCAC phenomenologically using low-Q 2 cross section data from 
the CHORUS experiment fl (see Sec. liVrl . 

III. NUCLEAR ASPECTS OF NEUTRINO SCATTERING 

In this Section we summarize the theoretical framework used in this paper to calculate the 
(anti)neutrino structure functions and cross sections for nuclear targets. A similar approach 



was recently applied in the studies of charged-lepton nuclear DIS 14J. In Sec. IIII Al we 
review calculations of nuclear structure functions in the impulse approximation with off- 
shell corrections, in Sec. IIII Bl we discuss nuclear pion correction, in Sec. IIII CI we deal with 
nuclear effects from coherent nuclear interactions of hadronic component of an intermediate 
boson, which are relevant at small x. 



1 To this end a covariant form of longitudinal vector = (<7 M + p^Q 2 /p ■ q) /Q is useful. 
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A. Incoherent scattering from bound nucleons 



In the region of large Bjorken x the nuclear DIS is usually treated in impulse approxi- 
mation, i.e. as an incoherent sum over bound protons and neutrons. In this approximation 
the nuclear structure functions can be written in terms of the convolution of the nuclear 
spectral function and off-shell nucleon structure functions. In the target rest frame we have 
(for a derivation and more details see 14, [32|, 33, 34]) 



F*(x,Q 2 ) 

F l( X iQ 2 ) 

l 2 F 2 A (x,Q 2 ) 
xF 3 A (x,Q 2 ) 



E jlMV T (e,p) (l + ^) (f t 

Y, j[Mv T (e,p) (i+^) (n 

"=p,n 

E fmr&p) (i+^ 

"=p,n 



Q 2 2 

Q 2 2 
6x /2 p^_ 



Q 1 



r 2 ■ 



Pz 
7M 



x'f;. 



(23a) 
(23b) 
(23c) 
(23d) 



where the integration is taken over the four- momentum of the bound nucleon p = (M+e, p) 
and [dp] = de d 3 p/(27r) 4 . The axis 2; is chosen such that q = (qo,0±,—\q\), p± is the 
transverse nucleon momentum, and 7 = \q\/qo- The energy and momentum distribution of 
bound nucleons is described by the proton (neutron) nuclear spectral function V p( - n \e,p) 
which is normalized to the proton (neutron) number in a nucleus. In the integrand Ff are 
the structure functions of bound proton (neutron), which depend on the Bjorken variable 
x' = Q 2 / (2p ■ q), momentum transfer square Q 2 and generally on the nucleon invariant mass 
squared p 2 = (M + e) 2 — p 2 . The relation between the Fx, F^ and F 2 of the off-shell nucleon 
is similar to (TUD, i.e.7 /2 F 2 = F T + F L with 7 ' 2 = 1 + Ax' 2 p 2 /Q 2 . Note that 7' = 7 for 
the vanishing momentum of the bound nucleon. It should be also noted that the transverse 
motion of the bound nucleon in the target causes the mixture of different structure functions 
in Eqs. fl23ajl and (l23bl) to order p 2 J(M\q\). 



1. Isoscalar and isovector contributions 

Since complex nuclei typically have different numbers of protons and neutrons, the nuclear 
structure functions generally have both the isoscalar and the isovector parts. In order 
to separate the isoscalar and isovector contributions in Eqs.( !23l . we write the isoscalar 
(P p+n = V p + V n ) and isovector (pp~ n = 7 ?p — V n ) nuclear spectral functions as 

V p+u = Av ^ ( 24a ) 

V p-n = (Z-N)V U (24b) 

where the reduced spectral functions Vq and V\ are normalized to unity. Using Eq. <124j) we 
explicitly write Eqs. (T2"3"|) in terms of isoscalar and isovector contributions 

K A /A = (F t N ) + I {F?- n ) 1 , (25) 
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where F t N = \{Ff + F"), Ff~ n = Ff — F™ for the structure function of type i and the 
parameter (3 = (Z — N)/A describes the excess of protons over neutrons in a nucleus. The 
quantities (F) and (F) 1 are the contracted notations of the integration in Eqs.(l23l) with 
reduced spectral functions Vo and Vi, respectively. The model of Vo and Vu which is used 
in numerical applications in this paper, is discussed in more detail in Ref. |l4j|. 

Equation fl25|) is generic and can be applied to any structure function. We now discuss 
in more detail the v ± v combinations of neutrino structure functions with definite C-parity 
(see Eq.(J7])). Let us first consider symmetric v + v combination. From Eq. fl25]) we have 

p(»+*)A /A = (iff +p )*^ + ((3/2) (iff +P)(p " n) ^ (26) 

for any type i of the structure function. 

In the absence of heavy quark contributions we have iff +l ^ p = _pf because of the 
isospin symmetry. For this reason the isovector term in Eq. (|26l) should vanish. However, we 
should remark that the isospin relations for structure functions for neutrino CC scattering 
are violated by the mixing of different quark generations and the c-quark mass effect, even in 
the presence of exact isospin symmetry on the PDF level. This effect results in the nonzero 
difference _pf qc sin 2 6q with 6q the Cabbibo mixing angle. Since the parameter (3 is 
small, this effect is suppressed in Eq. (l2"B"l) and, therefore, it is a good approximation to keep 
only the isoscalar term in Eq. ([26]) . 

Let us now discuss the v—v asymmetry in the nuclear structure functions. From Eq. (l25]) 
we have for the nuclear structure function of the type i 

Ft~ D)A /A = (^~ P)7V ) o + 03/2) (Ff -^-"^ . (27) 

The application of this equation to Fi and xF% requires somewhat more attention. We first 
consider F%~ v (a similar discussion also applies to Ft and Fl). This structure function is 
C-odd and dominated by the isovector quark distributions. In the absence of the Cabbibo 
mixing effect we have F^ = —iff ~ u ^ n . It follows from this relation that the first term 
in the right side of Eq. (]27p vanishes and the nuclear structure function is determined by the 
second (isovector) term. Nuclear effects in this case are illustrated in Fig. [TJ 

However, the isospin relations for structure functions are violated by the heavy quark 
effect and the Cabbibo mixing angle, as was discussed above. This effect generates a nonzero 
value of iff V ^ N oc sin 2 $c- Surprisingly, this effect should not be neglected in the analysis 
of nuclear corrections for v — v asymmeries even in the first approximation. This is because 
the isovector contribution in Eq. (l2"T|) is suppressed by the factor (3 and thus the relative 
effect of i<f is enhanced. This effect is further discussed in Sec. IIVt where we present 
the results of numerical calculation of nuclear structure functions. 

Let us discuss Eq. (I27|) in application to xF 3 . Note that xF^~ v is C even and includes 
the contribution from the light quarks, which determines the isovector nuclear correction 
in Eq. (1271) . and the s-quark contribution to the isoscalar part in Eq. (l2"Tj) . In fact, the 
difference xF^~ v ' N is driven by the s-quark distribution [recall that in the parton model 
xF^ v ' N = x(s + s)]. Because (3 is a small parameter, the relative contribution from the 
isoscalar term in Eq. (!27l) is enhanced and we conclude that the asymmetry xF^~ v in nuclei 
is dominated by the strange quark contribution. 
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FIG. 1: The isovector nuclear effect (F^ P)p )i/F^ P)p calculated for 207 Pb at Q 2 
and without off-shell correction using the isovector nuclear spectral function of Ref . [141] . The labels 
on the curves correspond to effects due to Fermi motion and nuclear binding (FMB) and off-shell 
correction (OS), which was calculated using Eq. (|30l) , 

2. Off-shell effects 



The structure functions of the bound proton and neutron can differ from those of the 
free proton and neutron. In the approach discussed, this effect is related to analytical 
continuation of the structure functions to the off-shell region and the dependence on the 
nucleon invariant mass p 2 . One can separate the two sources of p 2 dependence: (i) the 
target mass correction leading to the terms of the order p 2 /Q 2 and (ii) off-shell dependence 
of LT structure functions (or parton distributions). The off-shell dependence of the target 
mass effect is evaluated by the replacement M 2 — > p 2 in the expressions of Ref. [22| . In order 
to facilitate the discussion of the off-shell dependence of the LT structure functions, we note 
that characteristic momenta and energies of the bound nucleon are small compared to the 
nucleon mass, and the nucleon virtuality v = (p 2 — M 2 )/M 2 is a small parameter. We then 
expand the LT structure functions in this parameter and keep the terms to the first order 
in v, 

F 2 LT (x, Q 2 , p 2 ) = F 2 LT (x, Q 2 ) (1 + Sf 2 (x, Q 2 ) v) , (28) 
Sf 2 = <91nF 2 LT /<91np 2 , (29) 

where the first term on the right in Eq. (l28l is the structure function of the on- mass-shell 
nucleon and the derivative is evaluated at p 2 = M 2 . Similar expressions can be written for 
other structure functions. 

The off-shell correction 5/ 2 was studied phenomenologically in [3] by analyzing the data 
on the ratios of structure functions of different nuclei (EMC effect). As a result it was found 
that 5/2 is independent of Q 2 with a good accuracy and can be parametrized as 

5/2 = C N (x — xi)(x — x )(l + x — x) (30) 

with the parameters extracted from the fit Cn = 8.10 ± 0.30(stat.) ± 0.53(syst.), xq = 
0.448 ± 0.005(stat.) ± 0.007(syst.), and x x = 0.05. 
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Although the off-shell function (129!) generally depends on the structure function type, it 
was assumed in Ref. [3] that the relative off-shell effect is common for all types of structure 
functions (or partons) and described by the off-shell function (|30|) . It was tested that this 
assumption leads to fairly good cancellation between off-shell and nuclear shadowing effects 
in the normalization of the valence quark distribution. In this paper the off-shell function 
by Eq. (l30l) is applied in our calculations of nuclear (anti)neutrino structure functions. The 
interplay between off-shell (OS) and nuclear shadowing (NS) corrections in the DIS sum 
rules is further studied in Sec. IIIIEI and HHP} where we discuss the Gross-Llewellyn-Smith 
and the Adler sum rules for nuclear targets. 

B. Nuclear pion correction 

The lepton scattering off the meson fields in nuclei gives rise to mesonic contribution to 
the nuclear structure functions (for a review see, e.g., [H, [II]). We note that because of 
binding the nucleons do not carry all of the nuclear light-cone momentum and the meson 
contribution is important for balancing the missing momentum. At high Q 2 the mesonic 
correction to nuclear structure functions can approximately be taken into account using the 
convolution model and the kinematics of the Bjorken limit neglecting power corrections. 
Then the pion correction to F^ can be written as 

6 n Ff(x,Q 2 ) = Y,( &yh/A(y)FZ(x/y,Q 2 ), (31) 

where f n /A{y) is the distribution of nuclear pion excess 2 and F$ is the pion structure function 
and the sum is taken over different pion states n = 7r + ,7r~,7r°. Similar equations can be 
written for the pion correction to Ft and Fl with the same pion distribution function. We 
also note that the pion correction vanishes for F$. In the numerical analysis we also assume 
no off-shell effect in the virtual pion structure function, which would be the second order 
correction because (13T!) is a small correction by itself. 

We apply Eq. (|3~T|) to v ± v combinations of neutrino structure functions. It is convenient 
to separate the isoscalar and isovector contributions in the sum over pion states in Eq. (l3Tj) . 

We also apply the isospin relations for the CC neutrino-pion SFs (F^^ is isoscalar and 
independent of the pion charge state, and i^ - ^ 71 " = _p^ u ~ u ^ )_3 ^y e ^hen have for the 
nuclear pion correction 

^it + ^(*,Q 2 ) = f dyf° /A (y)Ff )n+ (x/y,Q 2 )i (32a) 

J X 

6 n Ft P)A (x,Q 2 ) = I dyfl /A (y)Ft P)7r \x/y,Q 2 ), (32b) 

J X 

where f®, A = f n + /A+fn-/A+fir° /a and f^ A = f n + /A—fw-/A are the isoscalar and the isovector 
pion distributions in a nucleus. In the applications discussed below, we neglect possible 



2 The contribution from the nucleon meson cloud should be subtracted since it is accounted for in the 
nucleon structure function. 

3 We neglect here the violation of the isospin relations because of the Cabbibo mixing angle (see Sec. lIII A"Tj) . 
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isovector contributions and use a model isoscalar pion distribution function, which obeys 
the constraints due to the nuclear light-cone momentum balance equation and equations of 
motion of the pion field in the nucleus (for more details see 14, 32|). 



C. Coherent nuclear effects (shadowing and ant ishad owing) 

At small x, at which the DIS longitudinal distance Lj = l/(Mx) is large compared to 
typical nuclear sizes, the nuclear DIS can be viewed as a two-step process in which the 
intermediate boson first fluctuates into a quark- ant iquark pair (or hadronic configuration) 
which then scatters off the nucleus. As an average lifetime of this fluctuation is large (of 
order Lr), it can undergo multiple interactions with bound nucleons (for a recent review see, 
e.g., [l2|). In this section we discuss nuclear effects in neutrino scattering associated with 
this mechanism and derive corrections for different SF types for v and v interactions. 

The structure functions at small x can be presented as a superposition of contributions 
from different hadronic states the virtual boson fluctuates to. For the helicity structure 
functions W and W±, as defined in Eq. (fTT|) . we have 

W h = '%2w h (i)<T h (i,s), (33) 

i 

where ah is the total cross section of scattering of the hadronic state i with the given 
helicity h — 0, ±1 off the target nucleon (or nucleus) with the center-of-mass energy s = 
Q 2 (l/x — 1) + M 2 , and Wh describes the weight of the given hadronic state in the wave 
function of the intermediate boson. 

We first consider the transverse structure functions with h — ±1 and approximate the 
sum over hadronic states in Eq. (1331) by a factorized form [lij 

W h (x,Q 2 ) = w h (x,Q 2 )a h (s) J (34) 

where ah is an effective cross section averaged over hadronic configurations and Wh is the 
remaining normalization factor for the given helicity h. At low Q 2 and small x the effective 
cross section can be approximated by average vector meson cross section in the vector meson 
dominance model (VMD) (see, e.g., [13]). As Q 2 increases, the averaging in (1341) involves 
the rising number of active hadronic states. 

In the longitudinal case (Wo) we explicitly separate the PCAC term, Eq. fl22|) . and assume 
that the rest, Fl, can be expressed by an equation similar to Eq. (1341) with an effective 
longitudinal cross section. 

In this paper we are concerned with the calculation of nuclear corrections to the struc- 
ture functions and cross sections. In the following we will assume that the normalization 
factors Wh and the pion decay constant /„- are not affected by nuclear effects. Then for the 
transverse structure functions with h = ±1 the relative nuclear correction is determined by 
the corresponding correction to the effective cross section, 

AT? sw £ - ^ rvrt 

h ZWZ + NWZ Za p h + Naf { } 

where SW^ is the nuclear structure function of helicity h with the incoherent term sub- 
tracted, 5Wh = — ZW% — NW£, and a similar definition is used for 5a^. For the 
longitudinal structure function the corresponding relation to the cross section will be dis- 
cussed in Sec JIII C4l 
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1. Multiple scattering corrections to helicity amplitudes 
In order to calculate nuclear modifications of effective cross sections 5a we apply the 



multiple scattering theory [35|, [36]. Let a' p ^ be the proton (neutron) effective scattering 
amplitude in forward direction. We will choose the normalization of the scattering amplitude 
such that the optical theorem reads Imo(s) = o~(s)/2 with a the total cross section. For the 
following discussion it is convenient to write the forward scattering amplitude as 

a=(i + a)a/2, (36) 

where a = Rea/Ima. The nuclear scattering amplitude a A can be written as 

a A = Za p + Na n + 5a A , (37) 

where the first two terms on the right side are the incoherent contributions from bound 
protons and neutrons. 4 The term Sa A is the multiple scattering correction (MS), which can 
be written as 

5a A = i j d 2 bdz 1 dz 2 e ikL{zi ~ Z2) a ■ p(l) a ■ p{2) e s ^ 2 ' a \ (38) 

J Zl<Z 2 

5(1, 2, a) = i / dz'a- p{b,z'), (39) 
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where a ■ p = a v pp + a n p n and p p (p n ) is the proton (neutron) nuclear density normalized 
to the proton (neutron) number and the integration is performed along the collision axis, 
which is chosen to be the z axis, and over the transverse positions of nucleons (impact 
parameter b). We also use a contracted notation of the nucleon position on the collision 
axis p(l) = p(b, zi ), etc. The exponential factor in Eq. (l38l accounts for multiple scattering 
effects (see, e.g., |37J) and for double scattering approximation e s should be replaced by 1. 

Amplitude (1371) describes the scattering of virtual hadron configuration of the interme- 
diate boson W* off a nucleus. In intermediate states in the MS series one considers the 
scattering of on-mass-shell hadronic states, and for this reason there appears a nonzero lon- 
gitudinal momentum transfer k^. The quantity L c = 1/ki is the longitudinal coherence 
length of this state. If m is the mass of intermediate state, then ki = Mx(l + m 2 /Q 2 ). 
Since the cross section in Eq. (I34p reffers to the effective intermediate hadronic state, we 
treat m 2 = m 2 ^ as a parameter. 

Note that Eq. (l38l was derived in the optical approximation assuming that the nuclear 
wave function factorizes into the product of the single particle wave functions. In this 
approximation possible effects of correlation between bound nucleons are neglected. We 
comment in this context that the correlations are relevant only if the coherence length L c 
is comparable to the short-range repulsive part of the nucleon- nucleon force, which is about 
0.5 Fm. However, this region is limited to relatively large x at which shadowing is a small 
correction (see discussion in Ref. [13]). The transverse momentum dependence of elastic 
scattering amplitudes was also neglected, that is justified by a small transverse size of the 
meson-nucleon amplitude compared to the nuclear radius. 



4 The incoherent term is discussed in detail in Sec lIII Al together with nuclear binding, Fermi motion and 
off-shell effects. Here we focus on multiple scattering corrections. 
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We now apply Eg. (158]) to describe nuclear effects generated by scattering of charged 
intermediate W + and W~ bosons. Let us first consider the contributions from the light (u 
and d) quarks. In this case the scattering amplitude obeys the requirements of the isospin 
symmetry and we can write the relations a(W ± p) = a{W^n) = a°± ^a 1 , where a and a 1 are 
the scattering amplitudes corresponding to the isoscalar and isovector nucleon configuration. 
Using these relations we have 

aiW^ -p={a°± ^a 1 ) p A , (40) 

where Pa = P P + Pn is the (isoscalar) nucleon density and (3 = (p p — p n )/ Pa is the relative 
proton-neutron density asymmetry. In order to facilitate the discussion of isovector effects, 
we assume similar shapes of the proton and neutrion densities. Then (3 = (Z — N)/A is 
independent of the nucleon coordinate for nuclear interior. We also benefit from the fact 
that (3 is a small parameter for heavy nuclei, and expand the multiple scattering series fl38|) 
in (3 and keep the terms to first order in (3. Then for the scattering of W + and W~ bosons 
we obtain 

5a A {W ± ) = T A {a°) ± lf3a% A (a°), (41) 

where 



T A (a) 



ia C 



In A, 



da 



2iaCn 



In A, 



a 2 C 



(42a) 
(42b) 



and the quantities C A 3 incorporate the multiple scattering effects and read as follows: 

C A {a) = J d 2 6d2 1 d2 2 e ifci(2l - 22) pA(l)pA(2)e 5(1 ' 2 ' a) , (43a) 
C A (a) = -llC A {a) = 

I d 2 &d^d^ 3 e lk ^-^p A (l)pA(2)p A (3) e s ^ a \ 

J Zi<Z2<Z3 

where the e s factor is given by Eq. (!39l) with a ■ p replaced by apA- The rate of multiple 
scattering interactions is controlled by the value of the mean free path of hadronic fluctuation 
in a nucleus If = (p^a)" 1 . If If is large compared with the nuclear radius (i.e. at low nucleon 
density or/and small effective cross section), then the e s factor can be neglected and the 
coefficients C A and C A determine the double (~ p 2 ) and the tripple (~ p 3 ) scattering terms 
in the MS series. If Zy is small enough then the e s factor should be taken into account. 

We find from Eq. (1411) that nuclear correction to the sum of W + and W~ scattering am- 
plitudes is determined by the isoscalar amplitude a , and the nuclear W + — W~ asymmetry 
is proportional to (3a l . Note that this asymmetry also depends on a through multiple scat- 
tering effects. The implications to (anti)neutrino structure functions are discussed in the 
next section. 



2. Dependence of nuclear effects on C -parity and isospin 



Up to now the helicity dependence of nuclear correction was implicit. In this section we 
apply the results of Sec lIIICTl to the combinations of helicity amplitudes which correspond 
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to the structure functions of interest. Recall that F T is given by the average (W + i + W_i)/2, 
and F 3 is determined by the asymmetry W+i — W-\. We will assume that the helicity is 
conserved in the MS series so that Eq. (l38|) generalizes to scattering amplitude with given 
helicity ah- 5 

The neutrino and antineutrino cross section correspond to the interaction of W + and 
W~, respectively It will be convenient to discuss the v + v average and v — v asymmetry 
since these combinations have definite C-parity. We use the notation a l T = k( a +i + a -i) f° r 
the average transverse amplitude with the isospin / = 0, 1 and a A = \{a ! +1 — aL\) for the 
corresponding asymmetry. Note that the amplitudes a° and a\ are C-even, while a A and 
o)p are C-odd. 6 

Let us first consider the v + v sum of amplitudes and use Eq. (|4ip for helicity h — +1 and 
h = — 1. In particular, we consider the helicity average (T) and asymmetry (A). Taking 
into account that \a° A \ <C |a° | in the low-x region, we expand T A (a\ l ) about a° to order 
(a° A ) 2 and obtain 

5a^ +P)A = 2T A (a° T ), (44a) 
Sa^ +P)A = 2a A T A (a° T ). (44b) 

Using these equations we compute ratios (135]) for effective cross sections. We recall that 
in the Born approximation a^ u+u ^ A = A2a° for either helicity state T or A and, therefore, 
the ratio Eq.flM} is 5Tl u+p = lm5a^ + ^ A /(2Alma°). Using Eqs.flHD and (ESD we obtain 

5TZ U T +P = 4 Re [{i + a° T ) 2 C A (a° T )] /(2A), (45a) 
5U A +P = Im [{t + a° A )T A (a° T )} /A, (45b) 

where a\ = 2 Im is the effective cross section corresponding to the transverse isoscalar 
amplitude and a{ = Re a[/ Im a[. If the real part of the amplitude is small, then the MS 
correction is negative because of destructive interference of forward scattering amplitudes 
on the upstream nucleons that causes shadowing of virtual hadron interactions. On the 
other hand, if the real part is large then the interference in the double scattering term is 
constructive, which would lead to an antishadowing effect. It should be remarked that 81Z A +U 
(the relative nuclear correction to the structure function F^ +u ) depends on the C-even cross 
section a^. The result is also affected by the interference of the real parts of the amplitudes 
in the C-even and C-odd channels. If we only keep the double scattering term in Eqs. (I45aj) 



and f!45bj) . we arrive at the following relation: 



51l A +i? /5n» T +p = 2 * aA< ? T , (46) 



a 



T 



where ay and qa refer to I = amplitudes. Factor 2 in Eq. (1461) is a generic enhancement of 
the double scattering correction for the cross section asymmetry which has a combinatorial 
origin (in the context of F% it was discussed in (l4l . 38|). Ratio fT46|) can be futher enhanced 



if the real parts of the C-odd and C-even amplitudes have different sign (this is indeed a 
realistic case as will be discussed below). 



5 However, the helicity-flip processes can occur because of spin effects, which are not considered here. 

6 The symmetry of the effective amplitude a} h becomes apparent if we recall the correspondence with the 
structure functions: ~ Ima^, F^~" - Ima^, - Ima^, F^~" - Ima^. 
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We now discuss the v — v asymmetry. We apply Eg. (PUT) to the amplitudes and a&. 
Similar to the v + v case, we expand T A (a° ±1 ) about the average amplitude a T in a series in 
a° A up to the first order term. We obtain 

Sat P)A = P W(4) + \a\a\T 2 A {a%)] , (47a) 

5at D)A = (3 [a\T*{<%) + a T a° A T 2 A (a T )] , (47b) 

where T 2 A (a) = ^^(a) [see Eqs.PSl)]. Note that the first equation in fl47]) is C-odd and the 
second C-even. The quadratic terms a 1 a° in (14 7p produce a correction to the leading term of 
order |a^/a^|. Taking into account the correspondence between effective amplitudes and the 
structure functions, we find that by order of magnitude |aA/ a rl * s the ra ti° of valence to sea 
quark distribution in the proton. This ratio is small in the region of small x (x <C 0.1), where 
coherent nuclear effects are most important. For this reason in the following discussion we 
keep only the first term in the right side of fj47|) . In this approximation the structure of 
Eqs.flJT]) is similar to that of Eq. fl44bl) . 

Using Eqs. (j4"Tl) we calculate the relative corrections to the cross sections. For this purpose 
we recall that in the Born approximation a (yV ~ v ^> A — (Z — A r )a 1 for either helicity state T or 
A. Then the ratio 5n v ~ D = Im5a^-^ A /[(Z - N) Im a 1 ] can be written as 

5K T J = Im [(i + a TA )T A (a° T )] /A, (48) 

where a T and a\ are the Re / Im ratios for the amplitudes a T and a A , respectively. 

Note also that in terms of (I, C) classification 5TZ T +D = 6n^ + \ 511 A +p = STZ^-\ 
51Z T ~ U = 5TZ^~\ and 51Z U A U = 57^ 1,+) . Figure [2] illustrates the x dependence of the ratios 
5TZ^ I,cr> for different isospin and C-parity states. 




0.0001 0.001 0.01 0.1 0.0001 0.001 0.01 0.1 

Bjorken x Bjorken x 

FIG. 2: The ratio 5TZ {I ' C) calculated for different isospin and C-parity scattering states for 207 Pb 
at Q 2 = 1 GeV 2 using the parameters of the scattering amplitudes from Sec. IIII C 51 The labels on 
the curves mark the values of the isospin / and C-parity, (/, C). 



3. Coherent nuclear correction to Ft and F% 

We now apply these results to neutrino structure functions. We recall that the discussion 
given above applies to the u- and <i-quark contribution and make use of the isospin sym- 
metry relations for the effective amplitude. In order to take into account the strange quark 
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contribution, we explicitly separate the ligh-quark (u and d) and the s-quark contributions 
to neutrino structure functions i = T, 3 

F v +V = F {u+u)(0) + F (u+u)(s)^ (49a) 
F „_ p = F (u-u)(l) + F {u-u){s)^ (4gb) 

where the superscripts indicate the contribution from the ligh quarks with total isospin 
or 1, or the s-quark contribution. Although this separation is quite general, we illustrate 
the terms in Eqs. (149j) in the parton model. For the transverse structure function we have 
F jy+v){o) _ a -^ li _|_^_|_^_|_^ j p(y v >w = x (^ ( i_ u j r y i _ d). The s-quark contribution is relevant 

for the C-even p T v+y ^ s ^ — x [s + s). The C-odd F T V u ^ s ^ is proportional to s — s and can be 
neglected. Similarly, for xF^ we have xF^ u+u ^ ^ = x(u+d—u—d), xF^ u ~ u ^^ = x(d—u—u+d) J 
and xF% u = x(s + s). The C-odd xF^"^ is similar to F T U v ^ and can be neglected. 

In general, the s-quark terms should be corrected differently from those of light quarks. 
Taking into account Eqs. (1491) we have for coherent nuclear correction to the structure func- 
tions Ft and xF 3 per one nucleon 

5F {u+-u)A = 5n (o,-) F (y+*M^ (50b) 
sf (u-p)a = ( 50c ) 

5Ft~ P)A = P5n {1 ' +) Ft~ pm + 5n^ +) Ft~ P){s \ (50d) 

where the structure functions in the right side should be taken for the free proton and the 
ratios 51Z are given by Eqs. ( 1451) and ( 1481) . Note that the s-quark contribution in Eq. (l50dl) 
is isoscalar and for this reason it is not suppressed by the factor (3. From Eqs. (150j) we derive 
coherent nuclear corrections to the (anti)neutrino structure functions. We define the ratio 
51Z\(A/ N) = 5F" A / 5F" N describing the nuclear correction of the neutrino structure function 
of type % in the units of the isoscalar nucleon structure function F" N = (F" p + F" n )/2 [and 
similar definitions for the antineutrino ratio 5TZ^(A/N)}. We have 



2F£ 2Ft 

F (u-D)(l) p{u-u){s) 



STZf\A/N) = 5TZ^ ± + (5K^ - ^<°.+>)£l_, (51a) 



snf \a/n) = srm ± ^a.+)£L__ ± {5 n^ - ^<oH)£l_ (5 ib) 

where the sign +(— ) corresponds to neutrino (antineutrino). In numerical applications 
in Sec. IIVI we will assume similar nuclear effects for 1 = light quarks and the s-quark, 
5K^ = 6K(°> + \ 



4- Coherent nuclear corrections to Fl and F2 

We now discuss coherent nuclear corrections to the longitudinal structure function. Nu- 
clear modification of the PCAC term is driven by the corresponding correction to the virtual 
pion cross section, and we have SF^ CAC / F^ CAC = Sa^A./ '(Act^n) = 5TZ n (A/N). The nuclear 
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correction to the pion cross section can be analyzed along the lines discussed in Sec. IIII C II 
We consider the isoscalar iiN scattering amplitude a° = (a w + p + a 7T - p )/2 and the isovector 
asymmetry a\ = a^+p — a^-p. Nuclear corrections to these amplitudes are given by equations 
similar to (145 al) and (14"8|) . and we have 

SK = ^ Re [ (< + alf #(a°)] /(2A), (52a) 
^ = Im[(z + ayT^(a°)]/A (52b) 

where <r° is the ttN cross section in the isoscalar state and a n = Kea n /lma n for the pion 
amplitude in the corresponding isospin state. In numerical applications we calculate 51Z n 



using 7rN scattering amplitude of 31J (see Appendix IAJ . Since the PCAC relation involves 
the cross section of the virtual pion with four-momentum q, the longitudinal momentum 
transfer in Eq.^gJ is k\ = Mx(l + ml/Q 2 ). 

The ratios 57t n for the cross sections of tt + and tt~ mesons can be readily obtained from 
Eqs.( l52|) . and we have 

5n w± = snl ± (ps-Rl - snl) {(J ^~ a ^\ (53) 

where the sign +(— ) corresponds to 7t + (tt^) scattering. 

The nuclear corrections to the structure function Fl, which include the vector-current 
contribution and non-PCAC terms of the axial current, can be treated similarly to the 
transverse structure function Ft- In particular, we assume that the relative correction to Fl 
is determined by the effective cross section in the longitudinal state and 571 l = 5F L A /Fl = 
5olI°l- Taking into account both contributions, we have for the relative nuclear correction 
to Fl 

5K l {A/N) = r PCAC 57i n {A/N) + (1 - r PCAC )57t L {A/N), (54) 

where r PCAC = 7 3 F| >CAC /Ff / should be computed for the isoscalar nucleon. 

Nuclear correction to F 2 can then be calculated in terms of the corresponding corrections 
for Fl and F? using relation (jlip . We have for the relative correction 

/ a / a r\ 571$ + 71571* 
5TZ 2 (A/N) = T i + R , (55) 

where R = Fl/Ft, and Fl and F? are the structure functions of the isoscalar nucleon. 



5. Effective scattering amplitude and cross section 

We now discuss the model of the effective scattering amplitude for different isospin and 
helicity states a T h which is used in our calculations. This amplitude determines the rate of 
nuclear multiple scattering corrections as discussed above. In Ref . [3] the C-even amplitude 

was studied phenomenologically using charged- lepton data on the ratios FA/F^ , and 
it was shown that a good description of data for a wide region of nuclear targets can be 
achieved using a constant and the effective cross section parametrized as 

a T = a /(l + Q 2 /Ql). (56) 
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The values of the parameters were fixed to a = 27 mb and ay = —0.2 in order to reproduce 
the low-Q 2 limit [37|. The scale parameter Qq, which describes transition to high Q 2 , was 
determined to be Q 2 = 1.43 ± 0.06(stat) ± 0.20(syst) GeV 2 . 

In this paper we test the hypothesis that the C-even effective amplitude a T also applies to 
the (anti) neutrino scattering and determine F^ +u . The other amplitudes (a° A , a T , a^J are 
not directly constrained by the CL DIS data. These quantities determine nuclear corrections 
to different combinations of neutrino and antineutrino structure functions as discussed in 
Sec. IIII CI It is important to note that only Re/Im ratios of these amplitudes are relevant 
for the calculation of nuclear corrections to the (anti)neutrino structure functions. In order 
to evaluate Re/Im ratios we apply the Regge pole approach, which proved to be very useful 
in analyses of high-energy hadron scattering (see, e.g., j39|), and approximate the scattering 
amplitudes by a single Regge pole with proper quantum numbers. We recall that Re/Im 
of the Regge pole amplitude in forward direction is determined by its intercept a#(0) and 
the signature (the latter corresponds to C parity of the scattering amplitude) Re/Im = 
— (C + cos naji(0))/ sm7ia R (0). The Regge poles of lo, p, and A 2 mesons have appropriate 
quantum numbers (see Table [J as known from hadron scattering phenomenology [39| these 
Regge trajectories have an intercept close to 0.5). In our calculations of (anti)neutrino 
nuclear scattering discussed below, we use the amplitude a T derived from the analysis of 



charged- lepton nuclear data 14|. The Re/Im ratio for the amplitudes and a T are fixed 



from the studies of the Adler and the Gross-Llewellyn-Smith sum rules for nuclear structure 



Amplitude 




Regge pole (I ) 


u(0-) p(l-) A 2 (l+) 


Re/Im 


1.15 (1) 1.35 (1) (-1) 



TABLE I: Regge poles contributing to the effective amplitude in different isospin and C parity 
states. Shown are the values of the Re/Im ratio extracted from analysis of nuclear corrections to 
the Adler and the GLS sum rules (see Sec. IIII Dl and HTl E|) and those calculated with the intercept 
or(0) = 0.5 common for all Regge poles (in parenthesis). 



The effective cross section gl for longitudinally polarized intermediate sates, which de- 
scribes nuclear corrections to Fl (see Sec. IIII C4|) . includes contributions from the vector and 
the axial currents. We separate the PCAC contribution to Fl, Eq.fl22l. and treat it explic- 
itly as described in Sec. Ill Dl and IIII C 41 The remaining part Fl is described by the vector 
and non-PCAC terms of the axial current. We assume that Fl is similar to Fl in CL scat- 
tering. In order to quantitatively evaluate <tl, we assume the relation o'l/o't = Fl/ Ft = R, 
which allows to calculate ol in terms of ot and R. This relation follows from (1341) if the 
normalization factors are independent of helicity. This can be justified in the VMD model 
at low Q 2 [l4| and we also apply this relation to the region of high Q 2 . For Q 2 > 1 GeV 2 for 
R we use the ratio Fl/Ft calculated using phenomenological PDFs and HT terms of (24j . 
In order to obtain R at Q 2 < 1 GeV 2 we extrapolate from large Q 2 , assuming R oc Q 2 at 
Q 2 -0. 

It should be also commented that the nuclear MS effect in the transition region of x from 
0.01 to 0.1 and Q 2 below 1 GeV 2 is affected by the value of the mass of effective intermediate 
state m e ff through Icl [see Eq. fl38|) ]. In our studies of the charged-lepton data we used the 
value m c fj = 0.9M which nicely fits the shadowing data. We also apply this parameter in 
the calculation of the vector-current contribution in neutrino scattering. Apparently, for 
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the axial current this parameter could be different because of different set of intermediate 
states. In order to account for this effect we use a simple ansatz of rescaling m c g with the 
ratio of a\- and p-meson masses, m eS m ai /m p . Figure [3] illustrates the dependence of the 
ratios STZ^°' + ^ and 51Z^^ on the parameter m c g. Note that this dependence is noticeable 
for the transition region of 0.01 < x < 0.1 and for low Q 2 . 
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FIG. 3: The ratios 5ft(°>+) and 8U^~^ calculated for 207 Pb at Q 2 = 1 GeV 2 using the scattering 
amplitudes from Sec. IIII C 51 and different values of the parameter m e jf. 



D. The Adler sum rule for nuclei 

The Adler sum rule relates the integrated difference of the isovector combination -F 2 1/-i/ 
to the isospin of the target: 

rM A /M 

S A = dx Ff(x, Q 2 )/(2x) = 2 I z , (57) 
Jo 

where the upper integration limit, the ratio of the nuclear and the nucleon masses Ma/M, 
is the kinematical maximum of the Bjorken variable x = Q 2 /(2Mqo) for the nuclear target, 
and I z is the projection of the target isospin vector on the quantization axis (z axis). For 
the proton = 1, and for the neutron S% — — 1. In the quark parton model the Adler sum 
is the difference between the number of valence u and d quarks of the target. The Adler sum 
rule survives the strong interaction effects because of CVC. However, in the derivation of 
the Adler sum rule the effects of nonconservation of the axial current as well as the Cabibbo 
mixing angle are neglected (for more detail see, e.g., (l^|). 

For an isoscalar nucleus the Adler sum rule is trivial since the integrand in Eq. (l57|) 
vanishes. For a generic nucleus of Z protons and N neutrons the Adler sum rule reads (we 
consider the nuclear structure functions per one nucleon) 

S£ = {Z- N)/A = (3. (58) 

We now discuss in turn the contributions to S\ from different nuclear effects. Let us first 
consider the Fermi motion and binding (FMB) corrections to F^~ v from Eq. (1271) . We will 
consider the kinematics of high Q 2 and neglect power Q~ 2 terms in nuclear convolution. By 
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performing the direct integration by x of Eq. (!27|) in the case of F 2 and using the isospin 
relations for the proton and neutron structure functions, we obtain 



Sa = (3 J^ (Ft v)p ) i = P + 5S° S , (59a) 
5S ™ =P J!^ F ^ U)P (^ Q 2 )$f2(x) («>! , (59b) 



where = (p 2 — M 2 ) 1 /M 2 is the nucleon virtuality averaged with the isovector spectral 
function V\. Note that the FMB correction cancels out in (I59ap in the impulse approxi- 
mation. The integral in Eq. (159bj) describes the variation of the Adler sum for the off-shell 



proton [the structure functions in ( I59bl) should be taken for the on-shell proton]. In the 
derivation of the OS correction to the Adler sum rule ( I59bl) we assume a universal off-shell 
function Sf 2 , i.e. common for proton and neutron and independent from the probe (neutrino 
or antineutrino). 

The nuclear pion correction to the Adler sum rule is computed using Eq. (132bl) . We have 

5 SI = 2 (n^+fA - n^- /A ) , (60) 

where n^/A = f dy f^/A^y) is the average nuclear pion excess of the given pion type. These 
quantities also determine the pion excess correction to the total nuclear charge. From the 
charge conservation we have n n +/A = n^-jA an d, therefore, 5S\ = 0. 

The nuclear shadowing effect in the Adler sum can be computed using the results of 
Sec. IIII CI We integrate the shadowing correction 5F 2 ~ U [see Eqs.( l50~]) ] and obtain 

ssT = pj Q 1 ^Ft u)p (x,Q 2 )sn^-\ (6i) 

where 51Z^~^ is given in terms of 5TZ^'~^ and STZ^'^ by Eq. (l55|) . 

It follows from Eqs. (1581) and (I59al) that the total nuclear correction to the Adler sum 
rule should vanish. This requirement provides an important constraint on the isovector part 
of nuclear structure functions. We verify that our approach is consistent with the Adler 
sum rule and explicitly calculate the OS and NS corrections as a function of Q 2 using the 
off-shell function 5f 2 by Eq. (l30p and effective cross section (1561) derived from analysis of 

nuclear DIS with charged leptons 14 1. It should be remarked that the quantity 571^' ^ is 



also determined by a^, the Re / Im ratio of the isospin 1 transverse effective amplitude [see 
Eq.( )48l) ]. This quantity is not constrained by charged-lepton data. We assume a\ to be 
energy- and Q 2 -independent. In order to constrain this parameter we apply the requirement 
of minimization of the total nuclear correction to the Adler sum rule £Sa = 8S® S + 5S^ S . In 
particular, we start from a\ = 1, as suggested by simple Regge arguments in Sec. IIII C~5| and 
try to adjust this parameter in order to minimize 5Sa for all Q 2 . We found that a\ = 1.35 
provides a good cancellation between the OS and NS corrections (see Fig. H]). 

We also remark that phenomenological cross section (I56p should be applied in the limited 
region of Q 2 . This is because this quantity was constrained by nuclear shadowing data in 



DIS which are limited to the region Q < 20 GeV [14j. In order to evaluate the effective 
cross section at higher values of Q 2 we take the requirement 5S\ = as equation on the cross 
section and solve it numerically. The result of calculation for iron is shown in Fig. HI Also 
shown is the effective cross section obtained from a similar calculation for the normalization 
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FIG. 4: Relative off-shell (OS), nuclear shadowing (NS) and total nuclear correction to the Adler 
sum rule calculated for iron (left panel). The right panel shows phenomenological cross section by 
Eq. (l56p with the error band and effective cross section computed by requiring exact cancellation 
between OS and NS effects in nuclear Adler sum rule (1) and in the normalization of nuclear 
valence quark distribution (2) (see text for more details). 

of the nuclear quark valence number (see Sec. lIIIEl) . Both approaches give consistent results. 
We also examined this approach for different nuclei and found similar results for effective 
cross section at high Q 2 . 

For convenience, we provide a parametrization which reasonably fits the numerical solu- 
tion to the effective cross section for Q 2 > 4 GeV 2 

af(Q 2 ) = 0.59 + 18.48 ln(g 2 /0.37)- L85 , (62) 

where the cross section is in mb and Q 2 in GeV 2 . In applications discussed in Sec. IIVI for 
the cross section, we use Eq.(|56|) for Q 2 < 15 GeV 2 and Eq.(lrJ2l for higher values of Q 2 . 



E. The Gross Llewellyn-Smith sum rule for nuclei 

The Gross-Llewellyn-Smith (GLS) sum is the integrated structure function F£ +v , 

i i-Ma/M 

Sgls= g / d*iT^(*,Q 2 ), (63) 

where we write the GLS sum per one nucleon for a generic nuclear target. In the quark 
parton model the GLS sum gives the number of valence quarks (baryon number) of the 
target 5q LS = 3 4(j. However, in QCD the direct relation between the baryon current and 



Sols only holds in the leading twist and leading order in a s . In contrast to the Adler sum 
rule, 5*gls depends on Q 2 and is affected by QCD radiative corrections, target mass, and 
the higher-twist effects. 

In this section we discuss the GLS integral for nuclear targets. We explicitly separate 
nuclear corrections to the GLS integral as Sq LS = Sq LS +5Sgls, where 5q LS refers to the GLS 
integral for the nucleon (proton). We first consider nuclear corrections to the GLS sum in the 
LT approximation. The correction due to nuclear binding and Fermi motion calcel out in this 
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approximation as follows from the direct integration of Eq. (123dl) (see also [34|). The nuclear 
pion correction to F^ +v also vanishes. However, both, the nuclear shadowing (NS) and the 
off-shell (OS) corrections to GLS, are generally present. We compute the OS correction 
to F 3 using Eqs. (1251) and (}2"§j) and assuming common off-shell function Sf 2 (x) = Sf 3 (x) 
(see discussion in Sec. IIII A 21) . The NS correction is given by Eq.( l50bl) . We write these 
corrections for high-Q 2 kinematics and neglect power terms in the structure functions 

*<S&s = \ 4lt P)p (^ Q 2 )Sf(x) « , (64a) 

^gEs = lJ o dxF^(x,Q 2 )5TZp'\x,Q% (64b) 

where (v) = (p 2 — M 2 ) /M 2 is the bound nucleon virtuality averaged with the isoscalar 
nuclear spectral function Vq. 

The subscript LT indicates that the structure functions are calculated in the leading twist 
approximation in QCD. As mentioned above, in LO approximation in as, l^iur^ reduces 
to the valence quark distribution u va \ + d va \ and the GLS integral gives the baryon number in 
the target. For this reason the total nuclear correction to the GLS integral should cancel out. 
The phenomenological OS and NS corrections driven by Eqs. ( 1301) and ( 1561) indeed largely 
cancel in the GLS integral, as was shown in Ref.[ll| (see Fig. 5 in Ref.jllj]). It should be 
remarked that phenomenological cross section effectively incorporates contributions from all 
twists since it is extracted from data. Futhermore, at practice the application of Eq. (l56|) 
should be limited to Q 2 < 20 GeV 2 (see also Sec. IIII Dj) . Higher twists are known to be 
important at low Q 2 and for this reason we should not expect the exact cancellation between 
(IMal) and fjMbl) at low Q 2 . 

At high Q 2 the LO QCD approximation to structure functions becomes more accurate 
and the total nuclear correction to the GLS sum rule should cancel out for normalization 
reason. We evaluate the effective cross section, which provides this cancellation, from the 
equation <5iSgls = 0. This equation is solved numerically using the LO approximation for 
the structure function xF% and phenomenological off-shell function by Eq. (l30|) . 



It should be remarked that the quantity SIZ^' ^ , which determines the nuclear shadowing 
correction to the GSL sum rule, also depends on a^, the Re / Im ratio of the transverse C- 
odd and isospin effective amplitude [see Eq. fl50bj) ]. which is not well known. We assume 



a a to be independent of energy and Q 2 . In order to fix its value we follow an iterative 
procedure. In particular, we start from a° A = 1, as suggested by simple Regge arguments in 
Sec. IIII C 5[ and calculate effective cross section a^(Q 2 ) from equation 55*gls — 0. Then the 
value of a° A is ajusted in order to match (J^(Q 2 ) at high Q 2 with the cross section calculated 
from the Adler sum rule in Sec. IIII Dl We found that for a^ = 1.15 the two solutions are 
in a reasonable agreement. The results are shown in Fig. [4] (right panel) together with the 
results of similar calculation using the Adler sum rule. 



IV. NUMERICAL APPLICATIONS 

In this section we apply the approach developed in this paper to calculate SFs and cross 
sections for the targets relevant to recent neutrino scattering experiments. In Sec. IIV Al we 
present results for the PCAC contribution to the neutrino SF, in Sec. IIV Bl we discuss the 
results for heavy nucleus/nucleon ratios of different structure functions, in Sec. IIV CI we deal 
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with the numerical analysis of DIS sum rules for nuclear targets, and in Sec. lIVDl we present 
our results for nuclear differential cross sections. 



A. Low-Q 2 limit 



Figure [5] illustrates the magnitude of the PCAC contribution to F 2 for the nucleon and a 
few nuclear targets calculated by Eq. (!22|) . In this calculation we use the Regge parameteriza- 
tion of irN forward scattering amplitude and the total pion-nucleon cross section of Ref.[3l| 
(see Appendix |AJ . Note that the values of F 2 for heavy nuclear targets are systematically 
smaller because of the nuclear shadowing effect for the pion cross section (see Sec. IIII CI) . 
In Table Hi] we list the values of F 2 corresponding to the limit Q 2 — > and x — > for the 
same targets as in Fig. [51 Our results are consistent with the value extracted by the CCFR 
experiment on an iron target 
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Target 


i(p+n) 12 C 56 Fe 207 Pb 


CCFR ( 56 Fe) 


F 2 


0.325 0.268 0.235 0.204 


0.210 ±0.02 



TABLE II: The value of F 2 in the limit of vanishing Q for neutrino interactions on different 

-5 



targets. The numbers are extracted at x 
for comparison. 



10' 



The determination from CCFR 



46J is also given 



The PCAC term determines the low-Q 2 limit of Fl, and Ft should vanish as Q 2 in the 
limit of Q 2 — > 0. In order to describe the structure functions in intermediate region, we 
apply a smooth interpolation between the high Q 2 regime, which is described in QCD in 
terms of LT and HT contributions, as discussed in Sec. Ill CI and the Q 2 — > predictions 
derived from CVC and PCAC arguments. We use Q 2 m = 1 GeV 2 as the matching point 
between high- and low-Q 2 regimes. 

It is instructive to compare the low-Q 2 behavior of R = Fl/Ft for charged-lepton and 
neutrino scattering. In both cases Ft oc Q 2 as Q 2 — * 0. However, if Fl cx Q for the 
electromagnetic current, for the weak current Fl — > _F£ CAC , and thus Fl does not vanish in 
the \ow-Q 2 limit. Then the behavior of R at Q 2 <C lGeV 2 is very different for charged-lepton 
and neutrino scattering. In order to illustrate this effect we calculate R as a function of Q 2 
for two different x (the x-bins of neutrino data [1, 0, 0]) for the isoscalar nucleon and a 
number of nuclei. The results are shown in Fig. [61 



B. Nuclear structure functions 



Figure [7| shows the result of calculation of the C-even structure function F 2 +u of 207 Pb. 
The resulting EMC effect resembles that of CL F 2 (see Ref. 14| for more details). At large x 
the nuclear correction is driven by FMB and OS effects. We recall that the off-shell effect, 
which modifies the structure functions of bound nucleon, is important to change the slope 
of the EMC ratio, bringing it close to data. At intermediate x values we observe some 
cancellation between the nuclear pion excess and the nuclear shadowing correction. The 
latter is the dominant effect for x < 0.05, providing a suppression of the nuclear structure 
function. 
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FIG. 5: The PCAC term of the neutrino structure function F2 (7-F£ CAC ) calculated for x = 10 -5 
(left plot) and x = 10 -2 (right plot) as a function of Q 2 for a few different targets (labels on the 
curves). A value M PCAC = 0.8 GeV is assumed (see Section ITV|). 




FIG. 6: Comparison of the ratio R = Fl/Ft calculated for the isoscalar nucleon for the charged 
lepton (dashed-dotted line) and neutrino (solid line) cases at fixed values of x. The left panel 
corresponds to x = 0.045, and the right panel to x = 0.125. Also shown are the results for different 
nuclear targets ( 12 C, 56 Fe and 207 Pb from top to bottom). A value M PCAC = 0.8 GeV is assumed 
(see Section HV] for more details). Determinations from SLAC [43], CCFR [48| and CHORUS 
are given for comparison. 
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The corresponding nuclear corrections for the difference F 2 ~ u are shown in Fig. [8] after 
normalizing them to (3 = (Z — N)/A. There are two distinct contributions to the differ- 
ence F 2 U ~ U for nuclear targets [see Eq. (j27]) ]. The pure isovector contribution [second term in 
Eq. (j27|) ] is shown in Fig. [TJ A comparison with Fig. [H] indicates that the dominant contribu- 
tion to the difference F 2 ~ u for nuclei is actually coming from the isoscalar part [first term 
in Eq.( l27l) ]. which is due to the Cabibbo mixing and the heavy quark production effects. 

It must be also noted that coherent nuclear effects are particularly pronounced in the C- 
odd channel, as can be seen from Fig. [8] (the shadowing effect at small x). The antishadowing 
effect at x ~ 0.05 is a combined effect due to nuclear smearing (FMB) and a constructive 
interference in the coherent nuclear correction due to the real part of the effective amplitude. 

Using the approach described in previous sections we calculate the ratios 1Z 2 = 
F vA i F v(p+n)/2 and n » = F vA/ F v(P+n)/2 for QC neutrino inactions for the Fe and Pb 

targets. The results are shown in Fig. [91 The bulk behavior of 1Z 2 is similar to that of 
1Z 2 i n CL scattering (see [3]) for both the x and Q 2 dependencies. Note that for heavy 
nuclei with an excess of neutrons over protons TZ 2 > 1Z 2 for x > 0.1. This is because the 
neutron excess correction is positive for the neutrino case, while this correction is negative 
for CL scattering. Indeed, this correction is given by the second term in Eq. fj25j) . Because 
F ^(p-n) > q and pv(p-n) < q ? we obtain the different sign of the neutron excess correction in 
these cases. It should be also remarked that this correction is negative for antineutrino SFs, 
similar to the CL case. As a result there is a cancellation for the neutron excess correction 
in the sum F 2 +u , as can be seen in Fig. [7J It is interesting to note that in the antineutrino 
case different nuclear effects cancel to a high degree in the region x ~ 0.1 (see also Fig. ITTT) . 

It is also instructive to compare nuclear effects for neutrino SFs F2 and xF 3 . One observes 
from Fig. [9] that at large Q 2 (= 10 GeV 2 ) the ratios 1Z 2 and 1Z 3 turn out to be similar 
(although not identical) in the entire kinematical region of x. We recall in this context 
that the effect of coherent nuclear interactions at small x is quite different for F 2 and 
xF 3 (see Sec. IHI Cp . Futhermore, the nuclear pion correction to xF 3 cancels out, and the 
antishadowing effect for F 2 and xF 3 is generated by different mechanisms (the combined 
effect of the real part of the C-odd effective amplitude and the off-shell correction for xF 3 
and the interplay between pion and off-shell effect for F 2 ). Nevertheless, at high Q 2 the total 
nuclear correction to xF 3 is very similar to that for F 2 . 

The differences between TZ 2 and TZ 3 are more pronounced at lower Q 2 , indicating dif- 
ferent Q 2 dependence of nuclear effects for these structure functions. The shadowing and 
antishadowing effects for F 3 become stronger in this region as illustrated in Fig. O 

Also one should remark that the target mass correction to the structure functions strongly 
affects the EMC ratio in the region of large x and low Q 2 ~ 1 GeV 2 . This effect causes a 
noticeable Q 2 dependence of the EMC ratio at x > 0.6. We also note that in Fig. [9] the 
curves for Q 2 = 1 GeV 2 are given for the region W > 1.25 GeV and for this reason the region 
x > 0.6 is not shown. 

C. Neutrino DIS sum rules 

As discussed in Sections IIIIDI and IIII E[ we observe a remarkable cancellation between 
the off-shell and the nuclear shadowing corrections. This cancellation, which becomes more 
accurate at high Q 2 where the leading twist dominates, can be attributed to underlying 
symmetry through the conservation of the isospin (Adler) and the nuclear valence quark 
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FIG. 7: Different nuclear effects on the ratio ±F^ +P)A / F^ +P)p calculated for 207 Pb at Q 2 = 
5 GeV 2 . The labels on the curves correspond to effects due to Fermi motion and nuclear binding 
(FMB), off-shell correction (OS), nuclear pion excess (PI) and coherent nuclear processes (NS). 




FIG. 8: Different nuclear effects on the ratio \F% ■ > /(/3Ff P)p ) calculated for 207 Pb at Q 2 = 
5 GeV 2 . The labels on the curves correspond to effects due to Fermi motion and nuclear binding 
(FMB), off-shell correction (OS) and coherent nuclear processes (NS). 
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FIG. 9: Our predictions for the ratios of CC neutrino nuclear structure functions normalized to one nucleon and those of the isoscalar 
nucleon (p + n)/2 (left panel for Fi and right panel for xF%). The ratios were calculated for 56 Fe (upper panels) and 207 Pb (lower panels) 
as a function of x at Q 2 = 1, 5, and 10 GeV . 
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number (GLS) in nuclear targets. Assuming exact cancellation we evaluate the effective 
cross-section describing the nuclear shadowing at high-Q 2 . Both the Adler and GLS sum 
rules can be used for this purpose, providing consistent results. In this paper we choose the 
Adler integral since the corresponding isospin relation is accurate to all orders. Furthermore, 
the evaluation of the Adler sum rule for nuclear targets involves the off-shell function Sf2, 
which was determined from charged-lepton DIS data. 

We can then evaluate the GLS integral Sgls(Q 2 ) for different nuclear targets. Results are 
shown in Fig. [KB The nucleon integral Sq LS (Q 2 ) is calculated using the NNLO coefficient 
functions and the NNLO PDFs of (24[. We observe that the nuclear correction 5 Sols 
decreases progressively by increasing Q 2 and ranges from about 2% at Q 2 = 2 GeV 2 to 
less than 0.3% at Q 2 = 20 GeV 2 . The results for different targets are similar and noticeable 
differences are present only at the lowest Q 2 values. We also note the general Q 2 dependence 



for Fe is in agreement with the CCFR measurement [411] • The values of Sq LS and 5Sgls 



are listed in Table HTT1 for a number of nuclear targets and fixed value of Q 2 = 10 GeV 2 . 




FIG. 10: The GLS integral for different nuclear targets as a function of Q 2 . The nucleon xF% is 
calculated in LT NNLO approximation with the target mass correction as described in Sec. Ill CI 
The bottom panel shows variations with respect to the result obtained for average isoscalar nucleon. 



Data points are the CCFR extraction of the GLS integral for iron 41 1. 



Power corrections due to target mass effect are included in the calculation. Apart from 
the target mass effect the dynamical power corrections should generally be present. The 
contribution of such terms to the GLS integral was evaluated in (42| using QCD sum rule 
approach, predicting a negative correction of —0.1 at Q 2 = 3 GeV 2 . We also comment 
that phenomelogical studies of HT terms in F3 using CCFR data were reported in [25[ 
with, however, high statistical uncertainties. In this context we note the current analysis of 



combined (anti)neutrino cross-section data [26| that would allow to greatly reduce statistical 
uncertainty in the neutrino HT terms. 
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Target 


Comment 


Oct 




P 


ttpp nvntnn 


2 72fi 


n 


D 


full calculation 


2.717 


-0.009 


12 C 


full calculation 


2.707 


-0.019 


56 Fe 


full calculation 


2.710 


-0.016 


207p b 


full calculation 


2.717 


-0.009 



TABLE III: The GLS sum rule calculated for different nuclei at Q 2 = 10 GeV 2 . The nucleon 



structure functions were calculated in LT NNLO approximation using the PDFs of Ref.[2J]. We 
use the lower cutoff x m ; n = 10 -5 in the calculation of the GLS integral. 

The nuclear effects somewhat modify the Q 2 evolution of the GLS integral. This is mainly 
due to the explicit Q 2 dependence of effective cross section describing the nuclear shadowing 
effect. As can be seen from Fig. [10] the correction is negative and of the order of 2% at low 
Q 2 . This must be compared with the variation related to the running a s , which is about 10% 
between 2 and 20 GeV 2 . Nuclear corrections should be then taken into account for precise 



extractions of a s from the GLS sum rule (see also [43]). The effect of nuclear corrections 
would indeed reduce the measured value of the strong coupling constant, by reducing the 
Q 2 slope in the GLS integral. 



D. Comparison with cross section data 

Using the results on nuclear structure functions we calculate the neutrino and antineu- 
trino inelastic inclusive differential cross-sections. Figure [TT] summarizes our results for the 
corresponding nuclear dependence. We note a cancellation between different nuclear effects 
for antineutrino at x ~ 0.1 -j- 0.2. As a result, the antineutrino cross section is not very sen- 
sitive to the target material in this region. Note also the cancellation of nuclear effects for 
x ~ 0.3 for the isoscalar target ( 12 C). Nuclear effect in this region is driven by the neutron 
excess. 

We can then compare our predictions with the available data on differential cross sections 
for different targets. Table HVl summarizes the most recent and precise measurements on 12 C 
(NOMAD), 56 Fe (NuTeV) and 207 Pb (CHORUS). 7 It must be noted that when comparing 
with data the predictions have to include electroweak (EW) corrections related to virtual 
loop diagrams, soft photon emission and hard photon emission. The effect of such processes is 
to shift the measured kinematic variables and hence to modify the differential cross sections. 



To this end, in this paper we use the recent one-loop calculations of Ref. [44j (see also 
Ref. (45|). We also implicitly assume that EW and nuclear corrections factorize so that they 
can be applied to the final nuclear structure functions. 8 

In Figs. [T21 to dH we report the comparisons between mesurements and our calculations 
for several values of the (anti)neutrino energy ranging from 20 GeV to 170 GeV. In this 



7 We use the cross-section data for all our studies and not structure function data provided by experiments. 
This procedure minimizes potential biases related to different model assumptions used in extractions of 
structure functions by individual experiments. 

8 We note that the magnitude of EW corrections in some specific kinematic regions can be comparable with 
nuclear corrections. 
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FIG. 11: The ratio of (anti) neutrino nuclear and the isoscalar nucleon differential cross sections calculated for 12 C, 56 Fe and Pb targets 
at E = 45 GeV. The electroweak correction is taken into account using Ref. 44 1. 
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FIG. 12: Comparison of our predictions (open symbols) with NOMAD data (full symbols) for neutrino differential cross-sections on 12 C 
E = 20 GeV (left plot) and E = 60 GeV (right plot). See text for details. 
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FIG. 13: Comparison of our predictions (open symbols) with CHORUS data (full symbols) for neutrino (circles) and antineutrino (triang! 
differential cross-sections on 207 Pb at E = 35 GeV (left plot) and E = 90 GeV (right plot). See text for details. 




FIG. 14: Comparison of our predictions (open symbols) with NuTeV data (full symbols) for neutrino (circles) and antineutrino (triang 
differential cross-sections on 56 Fe at E = 85 GeV (left plot) and E = 170 GeV (right plot). See text for details. 
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0.05H-0.95 


1423 
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56 Fe 
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35^340 
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0.10^0.80 
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25^-170 


0.020-^0.65 


0.10^0.80 


607 



TABLE IV: The list of neutrino and antineutrino data samples (after analysis cuts) corresponding 
to the cross-section measurements used in this paper. 



calculations we use both the nucleon parton distributions and the HT terms obtained from 
the fits to charged-lepton DIS data 24|. In the presented results the LT structure functions 
were evaluated in the NNLO approximation. In order to evaluate the HT contribution to 
neutrino structure functions F 2 and Ft we use the phenomenological charged-lepton HT 
terms rescaled according to the corresponding ratio of neutrino to the CL structure function 
in the LT approximation. No HT contribution is included for xF$. A more detailed study 
of (anti)neutrino data in the context of QCD fits and of the extraction of HT terms will 
be published elsewhere (26[. Note that the approach discussed was successfully applied in 
the analysis of charged-lepton nuclear data (the nuclear EMC effect) [13]. The comparison 
with neutrino data provides an independent verification of our model and also allows a 
compatibility check between (anti)neutrinos and charged leptons. 

In general, we observe good agreement between data and our predictions for all nuclei. 
We note that this also includes the lowest Q 2 data points, which in the case of CHORUS are 
at Q 2 ~ 0.25 GeV 2 . The Q 2 dependence in the low x bins is also reproduced, as can be seen 
from the y distributions and by comparing different (anti)neutrino energies. The existing 
cross-section data in the low x and low Q 2 region support the presence of the PCAC term in 
the longitudinal structure function. A value of the scale M PCAC = 0.8 GeV seems to provide 
the best agreement with data. 

After verifying the consistency of our cross-section model with (anti) neutrino data, it is 
interesting to use the model to examine the compatibility of data from different experiments 
and nuclear targets. At small and intermediate values of x the three experiments are in 
agreement for both the neutrino and antineutrino samples. At large values of x > 0.45 the 
NuTeV data seem to be systematically above calculations for most neutrino energies. On 
the other hand, the corresponding predictions for 12 C and 207 Pb at large x are in agreement 
with NOMAD and CHORUS data, respectively. 

The low-Q 2 and low- a; behaviour of the (anti) neutrino cross section is dominated by the 
divergence of the axial current and PCAC relation. However, as discussed in Section III D\ 
the scale controlling the PCAC contribution to the longitudinal structure function is not 
well known. The use of heavy targets introduces additional uncertainties. The MINERVA 
experiment recently proposed at Fermilab and currently under construction, has the 
possibility to address PCAC effect (as well as other nuclear effects) in heavy nuclei. To 
this end, the possibility to have cryogenic hydrogen and deuterium targets would greatly 
enhance the corresponding physics potential. 
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V. SUMMARY 

We discussed inelastic inclusive scattering of high-energy (anti)neutrino off nuclei and 
developed a quantitative model for nuclear structure functions using an approach which 
takes into account the QCD treatment of the nucleon structure functions and addresses the 
basic nuclear effects including nuclear shadowing, Fermi motion and nuclear binding, nuclear 
pions and off-shell corrections to bound nucleon structure functions. 

The presence of both the axial and the vector currents in neutrino interactions results in 
contributions with different C-parity to neutrino cross sections. The interference between 
the vector and the axial current is described by the structure function F3 which is not present 
in charged-lepton scattering. We discussed in detail how the nuclear effects depend on the 
structure function type (F 2 vs xF 3 ) and on the C-parity of structure functions (C-even u + u 
and C-odd v — v combinations). 

The axial-vector current plays a special role in neutrino scattering in the region of low 
Q 2 and low x. In this region the cross sections are dominated by contributions from the 
divergence of the axial current which is linked to the virtual pion cross section via the Adler 
relation. Using PCAC we examined the derivation of low Q 2 and low x limit for neutrino 
structure functions and discussed the scale in Q 2 at which the Adler relation can be applied. 
We studied this problem phenomenologically using low-Q 2 and low-re neutrino cross-section 
data. 

We examined the Adler and the Gross-Llewellyn-Smith sum rules for nuclear structure 
functions. A remarkable cancellation between nuclear shadowing and off-shell corrections 
was found, underlying the conservation of isospin and valence quark number in nuclei. This 
fact was used to constrain the effective scattering amplitude controlling the nuclear shad- 
owing for different C-parity and isospin states. 

We applied our results to calculate nuclear structure functions for the targets and 
kinematical conditions typical for recent neutrino experiments. Our predictions for the 
(antihaeutrino inelastic differential cross-sections agree well with the recent data on 12 C (H, 
56 Fe|2|, and 207 Pb Q. 

We conclude by commenting that the nuclear corrections prove to be particularly impor- 
tant for QCD phenomenology of v{p) data on heavy targets, where they significantly affect 
the fit results [26|. 
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APPENDIX A: PARAMETERIZATION OF THE PION NUCLEON 

SCATTERING AMPLITUDE 



The parametrization of the pion-nucleon forward scattering amplitude based on a Regge 



fit to ixN scattering data is given in Ref. [31 

1 
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In Eq. (IAip the amplitude aJr corresponds to the pion coupling to the isoscalar (isovector) 



nucleon configuration and the sign +(— ) corresponds to the 7r H 
neutron scattering amplitude is derived from isospin relations a n ± 7 
are 



(tv ) meson. The pion- 
= <Vf p - The parameters 



X 
Y 1 
Y 2 



12.08 ±0.29, 
26.2 ±0.74, 
0.560 ±0.017, 



e 

m 

V2 



0.0933 ± 0.0024, 
0.357 ±0.015, 
0.560 ±0.017. 



(A4) 
(A5) 
(A6) 



We also note that the forward scattering amplitude is normalized as Ima^ = ov/2 with 
er^ the total cross section, and the units of the parameters X and Y" 12 are such that a n is 
measured in mb if s is in GeV 2 . 
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